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APPENDIX: MODELS 1-10 


LIST OF NOMENCLATURE 


a 

P 

A 

P. 

C 

P 

C 

V 

F 

F 

g 

H 

k 

K 

K 

P 

L 

m(2, t) 

m 
M 
n 
n ; 

X 
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three hourly index of geomagnetic activity 
("equivalent planetary amplitude") 

daily index ("daily equivalent planetary amplitude") 
specific heat at constant pressure 
specific heat at constant volume 

flux of the solar 10.7 cm radiation (lo"^^ W/m^ c/s) 

10.7 cm flux averaged over 5 solar rotations (or 
5 months) 

acceleration of gravity of the earth (cm sec ^) 

scale height (pressure scale height) (km) 

Boltzmann constant (1.3805 x lo"^^ erg C°K]~^) 

thermal co^nductivity coefficient 

three-hourly index of geomagnetic activity 
("planetary' three-hour index") 

characteristic scale length 

coefficient in the density-solar flux relation 
(equ. (13)) 

mass of molecules or atoms (g) 
molecular weight 

-3 

total number density (cm ) 

number density of molecules or atoms (cm ) 

—2 

pressure (dynes cm ) 



-3 -1 

Q heat source (erg cm sec ) 

R gas constant (8.314 x 10*^ erg [°K] ^ Mol”^) 

-3 

p density (g cm ) 

T temperature (°K) 

t local time (hours) 

t conductive time (in formula (10) only) 

2 height (km) 

2 ^ height of lower boundary (120 km) 
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1. Introduction 


During the last years the structure of the upper at- 
mosphere above 200 km and its variations have been thoroughly 
investigated, in particular by studying the orbital decay of 
a large number of satellites with perigee heights in the 
range from 200 to 1500 km. The changes in the orbital 
periods of satellites due to atmospheric drag have indeed 
been proved to be the so far most informative tool to study 
the atmospheric structure above 200 km and its time-depen- 
dent behavior. 

The existence of the following variations in the at- 
mospheric densities has been established for the time inter- 
val from 1958 through 1963 , which covers nearly a full phase 
of decreasing solar activity within the eleven-year solar 
cycle. The variations are: 

1) a diurnal variation with a density maximum 
at 14:00 hours local time and a minimum at 
04:00 hours local time, 

2) a variation with an average period of 27 
days correlated with the solar decimieter 
radio flux in the 3 to 30 cm wavelength range, 
a semi-annual variation with a pronounced 
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minimum in June-July and a less conspicuous 
minimum in December- January , 

4) a long-term variation associated with the 
eleven-year solar cycle, correlated with 
the solar decimeter flux averaged over 
several solar rotational periods, 

5) a variation correlated with geomagnetic 
activity. 

It will not be possible here to acknowledge the 
numerous scientific papers v/hich revealed and confirmed 
the effects mentioned above. We therefore refer the 
reader interested in the details and history of these dis- 
coveries to recent reviews by Jacch j.a ( 1963 ^,b,^) and by 
Kinq-Hele (1963b ) . Therein one v/ill also find a discussion 
of how these effects are thought to be related to the 
physical processes in the thermosphere, in particular, to 
the heating of the atmospheric gases due to conversion of 
solar energy (either electromagnetic or corpuscular) into 
heat, and to the heat conduction (see also MacDonald 1963 , 
and Harris and Pries ter 1962a ) . All recent theoretical 
investigations are based on the generally accepted assumption 
that diffusive equilibrium prevails in the thermosphere and 


lower exosphere for the main constituents. xMass spectrometer 
measurements with rockets revealed that diffusive separation 
sets in slightly above the mesopause (sea the review by 
^ hamp . l . o n. ♦ usually the height of the turbopause is 

taken for the onset of diffusive equilibrium. The turbo- 
pause, defined by the decrease of vertical turbulence, has 
been observed by cloud-release experiments on rockets generally 
at a height of about 105 km. For the lighter constituents of 
the atmospheric gases, the escape from the earth can distort 
the diffusive equilibrium for these particular constituents 
in the thermosphere, when the exospheric temperatures are 


a jid__Pg±tex§illl (I,2.P.L) and Nicolet and Kockart-« 
( 15 ^) have shown that this affects essentially only the 
vertical distribution of atomic hydrogen. 

The first edition of the COSP.\R International Reference 
Atmosphere (CIRA 1961) was the first international model to 
give an account of the diurnal variation of the thermospheric 
structure by presenting atmospheric parameters for the maxi- 
mum - at 14:00 hours - and the minimum - at 04:00 hours 
besides a diurnal average. Since then, however, it has 
become apparent that with regard to the tremendous changes 
of the thermospheric and exospheric conditions during the 


solsr cycle » any atmospneric model wliicli gives only aveirage 

conditions for a certain time or certain level of solar ac- 

cannot fulfil the requirements for a reference model, 
which ought to provide an easy comparison between measure- 
ments of atmospheric parameters taken at different heights 
and different times. 

Recognizing the need for a reference model which ac- 
coun uS for the essential variations in the structure param- 
eters, the COSPAR Working Group IV decided, during its 
meeting in Warsaw in June 1963, to construct m.odels for 
the thermosphere and exosphere in the range from 120 to 
800 km for every two hours in local time. The variation 

within the solar cycle shoulo be accounred ror by selecting 

10 different levels or solar activity. Since so far only 
density data from satellite drag analyses are available 
in sufficient number for constructing height- and time- 
dependent models, it was decided that the models should be 
primarily based on the observed densities. We used the 
density data (or models) derived by Jacchia and Slov;ey (1982 
and ISM) , King-Kele ( 1963a ) , May ( 1962 ) , Marov (1963 ) , 

P.aet 2 Qld ( 1963 ) , Roemer (1963 ) , and the recent density de- 
^c^^iri^tions from Explorer 17 by Brva.n t ( 1964- ) , Charr.'Dj..on et al. 


n 963 ) , and Slowev ( 1964 ) . For heights in the range from 


120 to 200 km, density information is still very scarce 
(see, for excimple, Mikhnevich 1964 ) . The models attempt 
to reach a reasonable agreement between the densities de- 
rived with the inflatable-sphere method below 130 km (see 
Champion ' s review 1963 ) , and the requirement to represent 
the satellite drag data above 200 km. In the same way, 
the boundary values for the number densities of N^, O^, O 
and He were chosen as a compromise between recent mass 
spectrometer data ( Pokhunkov 1963 . Schaefer and Nichols 1963 , 
Nipr et al. 1964 . and Spencer et al. 1964 ) and the necessity 
of representing the satellite densities. No attempt was 
made to change the boundary values with the level of solar 
activity during the eleven— year cycle since no observational 
evidence is as yet available. Temperature measurements by 
Blamont (1964 ) between 100 and 150 km indicate that there 
was no appreciable change during the years 1960 to 1964. 
These observations cover about half of the decreasing phase 
of the solar cycle. They therefore might be considered as 
a justification for using fixed boundary conditions for the 
temperature and the number densities at a height of 120 km. _ 
But a further investigation of this problem, in particular ' 
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' with a series of rocket experiments during a complete solar 

i 

I cycle is highly desirable. As far as a diurnal variation 

i 

of the boundary conditions is concerned, a recent theoreti- 
I cal investigation by Harris and Priester ( 1964 ) revealed 

j that diurnal variations’ of either the temperature or the 

I 

j 

j number densities cannot be invoked for explaining the di- 

I urnal behavior of the upper thermosphere (above 200 km) . 

I 

On the other hand, Hall et a . l . (1 963 ) and Kallmann~Bi';i 1 . 

and Sibley (1964 ) derived larger diurnal variations in the 
number densities in the altitude range below 200 km than 
one would expect from models wherein diffusive equilibrium 
: prevails above a certain turbopause height and where this 

height remains the same during the entire day. Since, how- 
ever, the available evidence does not seem to be fully 

! 

conclusive and sufficient, it was decided to maintain all 
boundary values at 120 km constant throughout the day and 
the solar cycle. 
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2 . Basic Relations 


The physical procedure used to represent the afore- 
mentioned data in their dependence on height and time is 
based on the usual assumption that the follov/ing physical 
laws are applicable and sufficient to describe the behavior 
of the upper atmosphere in the range from 120 to 800 km. 

The hydrostatic lav/ provides a basic relation between 
the pressure gradient and the density; 


do 

dz 


p g 


( 1 ) 


where p is the total pressure, p the density, g the accelera 
tion of gravity and z the altitude above the surface of the 
Earth. Both pressure and density will be functions of al- 
titude as v/ell as time. Equation (1) describes a static 
equilibrium. Nevertheless it will be applied here to des- 
cribe a quasi-hydrostatic equilibrium under the assumption 
that a sufficient hydrostatic balance is reached quickly 
at any time during the diurnal changes of temperature and 
density in the thermosphere. 

It can be expected that the atmosphere behaves like 
a perfect gas according to the equation of state 
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p = n k T 


( 2 ) 


where n is the number density, k the Boltzmann constant 
(1.38 X 10 erg/°K) and T the absolute temperature. As 
the upper atmosphere is a mixture of different constituents, 

V 

we define the partial density of the ith constituent by 
\ 


P 


i 


= n. m. 

X X 


(3) 


where n. is its nvimber density and m. its mass. In the 

following we shall either use the mass m^ (in grams) or 

the corresponding molecular weight based on the value 

12 

12.0000 for the C isotope of carbon. - 

.The total number density used in (2) is simply the 
summation over all constituents 


n = ^ n^ . (4) 

i 

It is assumed that diffusive equilibrium, prevails in 
the entire range of altitudes from 120 to 800 km for all 
constituents with the exception of hydrogen. Combining 
equations (1) and (2) , one derives easily the following 
expressions for the partial pressures p^, partial densities 
and partial number densities n^ for the different 
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constituents 


(z, t) 


(2, t) 


= P, (z^.t) exp J* ^32] 


1 o 


2q kT(z,t) 


(5) 


T(Zo,t) 


2 m.g(2) 




pj^ (2, t) 


T(z^,t) 


P • (2 , t) — f-. P 
1 O ^ T(z, t) L j k1 


2 m.g(z) 


kT(z,t) 


dz 


] (7) 


Herein the lower boundary in either the height where ‘ 
diffusive eguilibrium begins (below 120 to aooording to the 
assumed herein) or the reference heiqr 
The mean molecular weight is defined by 


conditions 


ignt z^ = 120 km. 


? n (2, t) M. 
M(2,t) =2 L i 


Z n^ (2,t) 
i 


( 8 ) 


ana the scale height H 


H = - -i ^ ^ Tfz.t) 

V p dz y M(z , t) g 


g (2) 


( 9 ) 


where R is the gas constant (R = 8.314 x 10’ erg deg'^ Mol'h 
The main feature of the thermosphere is the rapid in- 
crease of the temperature above the mesopause (in the height 
range from 100 to 200 to) and the leveling off of the 



temperature distribution towards an isothermal layer in the 
upper thermosphere (z > 300 km) and exosphere. This property 
of the thermosphere is caused by the combined action of 
the heating from solar energy and. of the heat conduction. 

In particular the rapid decrease of the conduction time 
increasing height, which is essentially due to the 
increasing mean free path length for collisions in the ex- 
ponentially decreasing density distribution, and the time- 
dependence of the heating due to the rotation of the Earth 
determine the time variations in the thermosphere. The . 
conduction time t is defined by 


t 


P c 


2L 


K 


( 10 ) 


where is the specific heat at constant volume, K the 
thermal conductivity, and L a characteristic scale length 
comparable to the scale height. in the upper thermosphere 
(z > 300 km), the conduction time is so short that a flat 
temperature distribution with height is reached at any 
moment regardless of the large diurnal temperature varia- 
tions of several hundred degrees Kelvin. On the other 
hand, at 120 km, the conduction time is large compared with 
one day. Therefore the large temperature gradient is 
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laaini.ained and i*: cannot change considerably at this height 
within one day, 

D-ae to the positive temperature gradient in the thermo- 
sphere, one does not expect any appreciable turbuletit 
vertical convection. One has, however, to consider the 
transport .of heat due to mass rfiotion when the diurnal density 
bulge buil.ds up in the morning and shrinhs again after 14:00 
hours local time. Furthannore the horizontal pressure gra- 
dients around the density bulg;'. v/ill cause large-scale hori-- 
zontal motions. But since the aorisontal scale heights are 
tv/o to three orders of magnitude larger thcui the vertical 
scal.e heights, it is not expected that horizontal convection 
caused by the pressure gradients has a considerable influence 
on the diurnal ?Dshavior of the thermosphere. More detailed 
inves tigation on horisontal convection and conduction, how- 
ever, seems to be desirable. 

An expression for the time-dependant energy balance 
lias been given by Harris and Pries te r (1 962a , 1963 ) which 
includes heat conduction and the energy transport due to 
the diurnal expansion or contraction when the density bulge 
builds up in the morning and shrinks after 14:00 hours local 


1 








IT r hr 


^T 


( 11 ). 


The first term is the vertical heat transport due to thermal 
conduction, the second term determines the energy transport 
due to the aforementioned expansion-contraction. The third 
term accounts for the different physical processes which 
provide the heat sources by absorption of solar energy 
and heat losses (taken as negative heat sources) due to re-- j 

radiation (airglow in the infrared and optical wavelength j 

range). In equation (11), K(T) is the thermal conductivity 

and C the specific heat at constant pressure. Harris and ■ 

P _ .. ! 

Priester carried out a simultaneous integration of equation 
(11) and (6) on an IBM 7090 computer with heat sources pro- 
vided by the absorption of solar EUV radiation by atomic 
and molecular nitrogen, and a heat loss due to the reradia- 
tion of atomic oxygen in the inrrared (transition from the 

to the level) . The temperature and the number densities 

at the boundary ( 2 ^ = 120 km) were maintained constant. It 
was shown that under these conditions the theoretical diurnal 
density variation could not be brought into agreemenu with 
the observed variation. In particular, the uheoretj^cal 
density maximum occurred at about 17:00 hours locaj. time 
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Con- 


rather than at 14=00 hours where it was observed, 
seguently an additional heat source C'second heat source") 
was introduced in such a way that agreement between the 
observed and the calculated densities and their diurnal 
variations was obtained. l„e real nature of this second 
heat source is not yet understood, it is. particularly, an 
open question as to what fraction of this heat source is 
required to account for the amplifications introduced in 
the theory rather than accounting for other energy sources 
than electromagnetic radiation from the sun. A recent 
inemenr of the theory which accounts for lateral heat 
conduction in a preliminary way and also for the hea. pro- 
vided by the absorption of solar radiation in the Schumann- 
Runge wavelength range did not change the necessary require- 

ment for the empirical correction term in the 

L-erui in nhe heat conductio 

equat-ion (Hgrris and Pr-i eQ-f-.=iv -r, 

. rurtnermore the affect 

of thermal diffusion for the licht^r 

^ae consricuents (Ke and H) 

has been included. 

Hevertheless. the computer program used by Harris and 
Priester provides a sufficient number of free para.meters in 
order to reproduce the observed densities. The program is 
-onvenient tool to fit the observed densities with a 


time-dependent model v/hich is in agreement v/ith the general 
physical laws and which also conserves the mass in a vertical 
cylinder above the height of 120 km during the diurnal density 
variation. 

For these reasons, COSPAR Working Group IV decided, at 
its meeting in Warsaw in 1963, to make use of this- computer 
program to reproduce the observed densities and, in parti- 
cular, the time of the diurnal maximum which, according to 
observations in the height range from 350 to 800 km, re- 
mained constant at 14 hours during the years 1958 through 


1963 (Jacchia 1964 ^ . 
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^ — table s of t ha atmospheric pro'oarties in the ra nge 
^o m 120 to 800_krg for ten different levels of solar 
activi ty 


In order to provide convenient interpolation, the 
Working Group decided to present the tables of the hetero- 
sphere ror 10 different levels of solar activity. As the 
index for solar activity, the flux of the solar 10.7 cm radia 
tion v/as cnosen, since measursraents at this wavelength have 
been carried out continuously since 1947 by the National 
Research Council of Canada. Jacchia ( 1963 ') and Paetzold 
have shov/n that both the 27— day variation and the 
solar-cycle variation of the atmospheric density and tem- 
perature can be correlated with the solar 10.7 era flux F, 
but the dependence on F is different for the two effects. 

For the long-terra solar cycle variation, an average flux is 
generally being used. As Roaraer ( 1963 ) has shov/n, a flux 
averaged over 3 to 5 solar rotational periods is appropriate. 
Figure 22 and Table 1 shov/ the raonthly averages of the 10.7 


cm flux for the years 1953 through 1S63 . Table 2 gives the 
averages over intervals of 5 months for the years 1953 
through 1963. We shall designate these 5-monthly averages 
with F and shall use these values to represent the smiooth 
long-term variation of solar ?.ctivity throughout this paper. 
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Solar flux at 2800 Mc/s (10. 
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In order to establish the solar cycle variation, the 

I 

observed densities had to be reduced for the 27-day varia- 
tion and for the relation with geomagnetic activity. Fur- 
thermore, the semi-annual variation had to be taken into 
account. For the reduction .of the 27-day variation in the 

i 

densities, an empirical relation derived by Roemer Q963 ) i 

has been used. It depends on height and local time (see 
also Priester 1961 ) . For the "geomagnetic" effect, Jacchia's Q 963 \ 
formula for the change of the atmospheric temperature with . 
the geomagnetic index has been used to convert tempera- 
ture changes into density changes by means of a theoretical 
model. For the most recent findings concerning the geo- 
magnetically quiet periods see chapter 5 of this paper. 

With respect to the semi-annual variation, it was decided ' 

to use the averages for the months September through Novombar 
as standards (see chapter 5) . These averages usually coin- 
cide with the values for the March maximum. 

It then was attempted to derive the models in such a 
way that their densities represent the atmospheric conditions 
for the following 5-month average solar 10.7 cm fluxes; 

F = 275, 250, 225, 200, 175, 150, 125, 100, 75 and 65. For 
the latter (F = 65) , an extrapolation had to be made which 
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is a' little uncertain since the solar aiininvura had not been 
reached by the end of the tiiaa-iricer /al covvsred by the 
observational data. A siri-dlar uncarte,inty exists for the 
high values (F s 225). 

Since the relation batv;een the atmospheric structure 
and the average 10.7 era flux could only be established for 
a large portion of the decreasing phase of solar activity 
and since for the forthcc.u.ng increasing phase the relation 
may not be identical, v/e shall roh: label the models v/ith 
the ave rage 10.7 cm flux, buc rather give tire relcrtion 
betv/een the mooex numbers anc me F values in a socarate 
table (Table 3) . Occasiono.lly s’-alx s.iso denote the' ' 
models 1 to 3 as models for lo v soler activicy (? = 55 ...100), 
the models 4 to 6 as models for medxux; solar activity 

~ 125. . . 17 5) , the raocels 7 co 1 as :xC'lels for "nigh solar 


activity (F — 200... 250; , r.r..:: ...ooel 10 for extrei'aelv high 


solar activity. 

The tables present the av.r.c spheric 


two hours of local time. The moceis anplv co the 


rrs ror every 
ecu a CO real 


and temperate zones of che earth. he attempe h 
made to include the dependence of x:he acrospi 
on the latitude, since the amount of informa ■ 


ip O X. Jb. 


avai.i.aoh.3 


3.U pIT03Gn'C XS iiO V0 


c r-';’*""' > x* y 


Table 3 


Relation between the model numbers and 
the 5-month averages of the solar 10.7 era flux F-' 


Model No 

F 

T (MAX) 

T (MIN) j 

A 

P 

At' (A ) 
P 

1 

65 

1024 

1 

705 j 

2 

19 

2 

75 

1064 

731 

2 

19 

3 

100 

1184 

805 

3 

27 

4 

125 

1318 

885 

i 

4 

34 

5 

150 

1460 

979 

6 

48 

6 

175 

1585 

1053 

8 

59 

7 

200 

1729 

1159 

10 

69 

8 

225 

1845 

1229 

11 

73 

9 

250' 

1969 

1317 

12 

- 77 

10 

275 

2083 

1413 

13 

81 


*derived from'data obtained during the decreasing phase of 
the solar cycle. The unit of F is 10“^2 w/m2 c/s. Further- 
more the maximum and minimum of the diurnal temperature 
variation is listed for each model. Also the average A 
value for each model is given in the next to the last ^ 
column. The last column gives the correction AT^(A ) which 
has to be subtracted from the nighttime temperature^T (MIN) 
in order to reduce it to Ap = 0, since in deriving the models 
a correction formula was used which is in the meantime 
superseded by new findings (sea chapter 5) . 
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The heights are given in steps of 10 hn rrom 120 to 
300 km, and in steps of 20 km from 300 to 800 km. The 
quantities listed are; temperature (^K) , density (g cm ). 
pressure (dyne cm7^) , scale height (km), mean molecular 
weight, and the nuuiber densities for 0^, 0, He and A 

(cm”^). Number densities for hydrogen are given for 
heights above 500 km only, since belov/ this .level the con- 
tribution of hydrogen to the acraospheric densroy is com- 
pletely negligible- The numerical values for the number 
densities of hydrogen at 500 km were chosen in general 
agreement v;ith the values given by - ‘ — 

(1963 ) . 

In the Figures 1 to 13 the dependence of the tempera- 
ture, the density, the mraber densities of O^. O, He, 

A and H, the mean molecular weight and the scale height 

on either local time or heighu is given ror 3 

models; No. 2 representing low solar activity (F = 7o) , 

No. 5 for mean activity (F = 15w) a^id -.o. - -*-0- 

the depen- 


solar activity (F - 250) . In L.n- 


■ ses wher. 


dence on height is shov/n, four values ot ,ic.-al u ■ 04, 

10, 14 and 20 hours - have been selecred as parameters. 
Figs. IS and 20 show, for every two neurs in -ocal 
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time, the relation between the exospheric temperatures 
and the 5-monthly averages F of the solar 10.7 cm flux 
which represent the long-term variation of solar activity. 
In Fig. 21 an example for the agreement between a few 
observational data and the corresponding models is 
graphically presented. Due to the complicated dependence 
of the atmospheric parameters on the different effects 
v/hich act in the upper atraosphere, it was not possible 
to show how the models represent the huge number of 
measurements, in particular, the numerous density values 
derived from satellite drag data. 


4 . The accuracy of the data 


In the following part, an atterapt will be made to 
evaluate some estimates about the accuracies of the listed 
quantities. If densities are derived from precisely ob- 
served changes in the orbital period of a satellite by 
accounting for how the' orbit inrersects the diurnal bulge 
("integration of the drag force along the orbit by means 
of an auxiliary model with diurnal bulge"), one might 
s^^^ly assume that the essential unc^ertainty lies then 

only in the product of drag coefficient C , area to mass 

D 

^3.tio A/m and in the accuracy or rhe auxiliary densi—y 
model being used. The precision of the orbital data, in 
particular of tne period, might be judged frora the correla- 
tion between the period variations and the solar 10.7 cm 

flux and the geomagnetic flux A or a , since it has been a 

P P 


lisned that the most precisely reduced drac 


:ta shew the 


best correlations v;ith the solar and geomagnetic indices. 
For spherical satellites, the uncertainty of the crag co- 
efficient is likely to be less than ± 10 percent. ?or non- 
spherical satellites, the uncertaimey or the abovementioned 
product Cj^*A/m might amount to ± 25 percent or more, i.n 
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particular if tV.a tunibling mode is unknown. A e-rxoiv3 
limi-cation on the accuracy car. exist if dsnsitxas eira ua- 
rived from oddly shaped satellites vjith perigee he3o-hr.s 
above 400 kin, since for those satellites the efracc of 
solar radiation pressure cannot be properly taken .U'^o 
account {Jfcchi a lif^.) . On the other hand, since -the mod 
densities are derived from several satellites, tn-a uiicei.- 
tainty of; the model will, in general, be smaller than 25 
percent. From P-oemer ' s cliEgram (I.9b2, I'xg. 5.1/) 
relation bstVv-esn average r:olar fluxes and model nusv^jers 
of the Harris-Priester (1962) models derived from seven 
satellites at different levels of solar activity, one 
learns that all carefully reduced data are scat c.. 
around a linear relation, in a vnr/ which correspo-nds to a 
maximum deviation of ± 25 percent in uhe .in .. j. vie-, 
sities (taken at a height of 500 ?<m) • Cor re 5ponoin-.j , 
the mean error can be estimated as i 10 percent,. 

The accuracy of the number aciisd. c- •> 

constituents V.^, snd 0 can onl^' b= ro>.;;jU.-V 

from the disagreement of the different mass speccrorac te.-. 
results and also frora the determination baa.ed o-n the at- 
tenuation of solar SUV radiation in the lower tntrvvw-pa; 



A comparison shows that a systeraatic uncertainty of a 
factor of two might still exist in some of the measured 
number densities. Since trie temperatures derived from 
the density distribution depend essentially on the nuraber 
densities of atomic oxygen, the systematic error of the 
absolute values of the exospheric temperatures might 


amount to ± 10 percent, although an error in excess of 5 
percent seems not very likely. This follows ilrorn an in- 
vestigation of the possibility of matching the observed 
densities using consioerably different boundary values 
for the main constituents ( Stein and Walk er 19S4) . The 
uncertaini_y of me temperature at 120 kra might be as large 
as ± 25 percent. 

it is likely that tiie teruperauure , the density, and 
consequently the number densities at 120 km undergo 
systematic changes within the 11-year solar cycle. These 
variations might even exceed the limits of = 25 oercent 
the temperature and uhe ractor or uwo for the density 
and number densities. It is urgently desirable to measure 
the atmospheric parameters in the lower thermosphere, in 
particular at the boundary height of 120 km, for an entire 
solar cvcle. 
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Another uncertainty exists in the choice for the 
height of the so-called "diffusion level" ("turbopause") . 
For calculation purposes it is generally assumed that 
mixing prevails in the atmosphere up to a certain height 
and that diffusive equilibrium exists above this height. 
This can certainly be regarded as an approximation only. 
Within this scheme it is essential to establish the proper 
height for this level. In the models presented here the 
level was chosen to be 115 km above the earth's surface 
for the heavier constituents (N^/ 0^# 0 and A) as a com- 
promise between direct density measurements and the mass 
spectrometer data (see the references in chaptex 1) . 

The uncertainty of the number densities of helium 
depends merely on the uncertainty of the altitude where 
diffusive equilibrium is reached for this constituent; 
furthermore it depends on the largely unknown thermal 
diffusion coefficient. A relatively low height is assumed 
for the onset of diffusion (IGO km) for helium. Thermal 
diffusion has been taken into account. The reason for 
choosing the rather low diffusion level was to provide a 
sufficient number of helium atoms at heights above 600 km 
for low levels of solar activity in order to account for 




the observed curvature in the density curves (see Kinc;-Hele 
Harris and Priester 1963 ) . It was concluded that no 
other constituents in an isothermal exosphere could provide 
the observed curvature. It could only be explained by a 
transition from the atomic oxygen layer to a layer v;here 
helium is the dominant constituent. All conclusions., 
therefore, on the helium number densities depend strongly 
on the validity of tlie observed curvature in the density 
distribution. 


Still very problematic is the accuracy of the number 
densities of atomic hydrogen, since there are no direct 
measurements available yet. Furthermore the combined 
action of escape and diffusion determine the total a.mount 
of hydrogen in the hetarosphera and the hydrogen distribu- 
tion with height. Kocbarts e nd Ficclet ( 1963 ) furnished 

detailed calculations on the hydrogen distribution based 

7 -3 

on an assumed boundary value of 10 era at lOQ hm. Since 
the escape strongly increases v.dth the exospheric uam- 
psrature, we find less hydrogen atoms at a height of 500 
km during periods of high solar acitivity than at the 
solar minimum. For all these reasons the number densities 


of hydrogen are highly uncertain 


meralv be 


- tn - 


and should 




regarded as an example to demonstrate their behavior during 
the solar cycle. Up to heights of 800 km the contribution 
of hydrogen to the total density is negligible, even for 
very low solar activity, when the diurnal average tempera- 
ture is about 860 Even for this low temperature, the 

contribution of atomic hydrogen to the total density is 
only a few percent. 

In the following table we summarize the estimated total 
ranges of uncertainty due to systematic errors for the 
main physical parameters. 

Total density : : (long term behavior for March or 

September-Deceniber averages (semi- 
annual effectO after short term 
variations (27 -day variation, geo- 
magnetic effect) have been reduced 
■ ■ ./ -for) ' 


altitude 

km 

■error , (percent) 
high solar activity 

error (percent) 
low solar activity 

120 

± 50 ■ 

± 50 

200 

. +40, - 10 

±25 

300 

+40,-10 

±25 

400 

±25 

± 25 

600 

± 25 

± 25 

800 

o 

± 25 
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Temper Ht-.nr/a 


altitude 

120 

200 

300 

400 

600 

800 


e rror (percent^ 
i 30 
i 20 
db 15 
i 10 
i 10 
± 10 


lJliab££-i3sr.si^ The range of uncertainty 

is represented by two factors with which the given numeri- 
cal values ought to be multiplied. The estimated factors 
a.ie for o^, O, He and A: 

0.5 3,nd 1 , 5 . 

As far as the accuracy of the hydrogen number densities 
IS concerned, an estimate is difficult to obtain, since 
there are no measurements yet for the densities of neutral 
hydrogen. Furthermore the escape of hydrogen is sansitivel 
dependent on the exospheric temperature. Therefore any 
systematic errors in the temperatures will change the ac- 
cumulation of hydrogen during the decreasing phase of the 
solar cycle. Thus the given number densities of hydrcc on 
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should be merely considered as a guide for the order of 
magnitude. This statement is in concordance with Kocharts 
(156.3.) , One can see from the new tables that 
even for very low solar activity with a diurnal average 
temperature in the exosphere of 870 °K the contribution 
of hydrogen to the density is negligible up to a height of 
800 km. 

The error values given above define the estimated 
total range; of uncertainly due to possible systematic 
errors for the long~term variation of the atmospheric 
parameters with the eleven-year solar cycle. If one, 
however, wants to compare’ individual measurements of at- 
mospheric parameters with the corresponding model values, 
one has to allow for the short-term variations described 
under items 2, 3 and 5 in the introduction to this section 
or in more detail in Jacchia 's (1963, 1964 ^ review. In 
the following chapter we shall describe the way to account 
for these effects. 
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5. 


W l t.^ individual measnrpin»nl-g ^ tlrosi-.lia,-.-„ 


For the comparison of model data with individual 
measurements of atmospheric parameters one has to account 
for the 27-day variation, the semi-annual variati^on and 
the geomagnetic activity effect, since the models provide 
only the smooth long-term variation of the upper atmosphere 
in its dependence on the time within the eleven-year solar 
cycle and the dependence on local time within the diurnal 
variataon. The long-term variation of. the atmosphere is 
related to the solar decimetric flux when averaged over 5 
solar rotational periods. For convenience we might use the 
5-monthly averages of the solar 10.7 cm flux F as given in 
Table 2. Therefore the atmospheric parameters for each 
individual value of F can be found by interpolating between 
the corresponding models according to Table 3. The exo- 
spheric temperatures in their dependence on F are also 

conveniently read off from Fig. ig and 20 for every even 
hour in local time. 

Relatively simple and sufficiently accurate formulae 
have been derived for the short-period variations of the 
exospheric temperature from the statistical analysis of 
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satellite drag measurements. The formulae given belov/ are 
based on the v;ork by Jacglii a, and S 1 ov;ey ( 1962 . _19^, 1 964 ) 
(see Jassllis. 2^^) • £aet^)2i3. Q.J2M.) , R oemer (1963) , May: 
(25_^) and Nswtoii_je_^a^^ (22J&4) . it seems most convenient 
to follow essentially the procedure proposed by Jacchla 
(2Si?.4) / • that is, to formulae as aivef) 

■feel.Q W„ t -3ie_ji2.^h^tti,TRg_jaii)j^u^ then to find 

t>Le_othexj^tffip^Dheric, ^ 

CQXSg.g£2a^Qg_^Q_J^.a_^a2^UaJfegJ,_jiigMt^^ temp e rature 


3i.) T he 2 7 -day variation 

-The difference betv/aen the exospheric temperature for 
a given value F of the solar 10.7 cm flux and the model 
temperature corresponding to F, the 5-raonthly 10.7 cm flux 
average, is, for the nighttime minimum (04 hours local 
time) , given by 


AT - 1.9 (F - F) C°K] ( Jacchia 1964 ) .(12) 

The delay time for the response of the atmosphere to 
variations in the heat source as represented by the values 
F of the 10.7 cm radiation has been found to lie between 
0.5 i 0.3 and 2.1 days (for details see Jacchla * s review 
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1364 and Bgnrdeaq J^) . it is therefore advisable 

to use the measured 10.7 cm flux of the preceding day, if 
one wants to calculate the teinparature or density for a 
particular day. 

If one is only interested in the densities one might 
use the relation 


^ P miz,t) (log F — log P) 

where m(z,t) is an empirical factor derived 
statistics of the 27-day variation of the densities by 
(1963) . The following values might be used 


(13) 



Factor m(z , t) i r 

A 

Ponnula 

_il3) 


Loc. time t 
(hours) 

! 

1 200 

1 

1 3 00 

c (km) 
! 

1 400 

i 

i 

1 500 

1 ^ 600 km 

0 

1 

CO 

0.3 

1 

1 

1.0 

i 

i 1.7 

i 

2.2 

1 

1 2.7 

8-12 ! 

' 

0.3 

. 85 

t 

1 X • 4c 

1.7 

! 2.0 

12-20 ’ 

1 

0.3 1 

1 

.7 

1.1 

i -1- . ^ j 

! 1 

f 

1.3 

20-24 j 

■ 

O 

.85 

1.4 

i 

1 

1.7 1 

i 

1 

2.0 


Tne accuracy of the coefficients in formulae (12) and 
exactly knov;n yet, since they might 


(13) is not too 



themselves depend on. the eleven-year solar cycle,. Fur the 
more, the relation between the 10.7 cm flux F and the re- 
lated heat source of the atmosphere miaht not be 
and might not remain exactly the same for a v;hol.s solar 
cycle. The heat source responsible for the 27-day varia- 
tion is genercilly believed to be the absorption of solar 
extreme ultraviolet radiation, raore or less exclusively. 
Ihe ObO-1 satellite rev-aaled a clos-e proportionality 
between the averaged EUV flux and the solar 10.7 c-.i flux 
for the spring of K- 52 (M eu riart et al. 1964 ) , From there 
considerations one can expect the error in the temperature 


correction AT to ba 1 

le tween i 

10 

and i 

: 20 

perr 

5nt. A 

rough estimate might 

a.lso be o 

■j,. J,. .. 

ill ad 

for 

each 

individual 

date from the change 

in the da 

ily 

f lime 

val 

■ues 

F on the 

preceding days, since 

it canno 

t h, 

a firi 


daci. 

and yet 

whether the response 

time of ti 

he i 

atm.os; 

pher 

e 13 

closer 

to two days or only a 

fev; hour; 

3 • 

Thus 

uiie 

dif . 

cerenca 

between the AT values 

obtained 

for the 

t\vO 

prac 

-c'ding d-ays 

reveals an individual 

estimate 

of 

the a ecu: 

racy 

for each 


case^ 


aeoTiacnet-j c activ;! tv effect 

In two recent analyses of the variations in tlie upper 
atmosphere and their relations to the indices of the geo- 
magnetic disturbances (vJ..gccni a and S1,owey 1964 , Newton 

I it was discovered tnat Ciuring quiet intervals 

(A < 20, K < 4) the hea 
P P 


between the additional hea 

the index a . Jacchia cri’ 
P 


,ng or 

tne ther 

ations 

of geora 

This 1 

OG to an 

ing of 

the upp 

s , i n 

OiiG up da’ 


itrnosphere and 


his review on atiTiospheric suracture (iJiM.) / the foll ,>wihg 
expression loi tne increase 'Che exospheric tair^parature 
with the three •‘"Hou rly inoex cr geoeicigne cic activity a 


AT 


1?C a 




exp ( — 0. u3 a. ) j 
P 


(14) 


Prom a discussion of the density gauge maasurerr.ents on 

Explorer XVII ( Newton et al. IW) and the drag data 

(SlQv^ey 19,64) , one obtains a somewhat sirapler analytical 

expression for the range from a = 2 to a = 150 The 

P P 

change of the exospheric temparatura is represented by 


AT = 125° 


aog 


r> 


(15) 
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The zero-point of AT in this relation is chosen for = 2 

This yields the factor 1/2 in formula (15). The scatter 

of the individual data around these analytical expressions 

is of the order of i 50 °K. 

The following table gives a comparison of the two 

formulae (14) and (15). One sees that in a wide range 

the agreement between the values if far within the ecatter 

of the individual data. As far as the range a < 2 is 

P 

concerned, there is no conclusive information yet on how 

the extrapolation to a “0 should be made and what kind 

P 

of analytical expression represents the real behavior best 

Table 5 


Comparison of the increase of the exospheric 
temperature as a function of ap according to 
formulae (14) and (15) 


a 

P 

equ. (14) 
AT 
C°K] 

equ. (15) 

at.: :" 

:■ C°K] 

a 

P 

equ. (14) 
at 
C°K] 

equ. (15) 
AT 
C°K] 

0 

0 ' 


30 

144 

147 

1 

11 

' ■ 

40 

160 

163 

2 

21 

0 

50 

173 

175 

3 

30 

22 

70 

195 

193 

4 

38 

38 

100 

225 

212 

5 

46 

50 

150 

275 

234 

7 

61 

68 

200 

325 

(250) 

10 

79 

87 

250 

375 

(264) 

20 

120 

125 
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The geomagnetic index used so far was the three-hourly 

planetary index a . Often it is sufficient to use the 

daily average S a^. One can easily see that the 

formulae (14) and (15) can also be used for A if the a 

P p 

values do non deviate too much from their average value 
during the 24-hour period, that is, if the expression for 
an average of the logaritnm log A = 'T L log a yields 
essentially an identical value for A as the linear average 
A = ^ S a . 

p 8 p 

How shall the correction for the exospheric tempera- 
ture be applied to the n;.odal te’u;peratures? Since the im- 
portance of the additional heating of the tharmosohere 
during small changes of a^ ac geonagnetically quiet 
periods v/as only recently discovered after the tables 
ior the 10 differenr. levels or solar activity 'v'.'ere already 
pi-oduced, it was nor possio_e in tnis edition to reduce 

all models to exactly the same value of a , for exeumole, 

o 

to = 2 or even a^ = 0. An e:ca.r.;'.ina,tion of the basic 
satellite drag data from which the models are derived 
showed that the models for high solar activity correspond 
to an average A level of about 10, whereas the models 
for low activity correspond to an averace A of abnu-f* 

D 

2» There is a graoual increase in the 
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A level of the models from model 1 to 10. The values are 

P 

given in the second to the last column of Table 3 in Chapter 

3. It must be left to a future improvement of the atmospheric 

tables to provide models all reduced to the same A level 

P 

using the improved relation given in formula (14) . In de- 
riving the exospheric temperature for a certain date repre~ 
sen ted by F, F and A , one has to allow for this shortcoming 
of the tables when correcting the exospheric temperature for 
the proper A^ value. After one has derived the exospheric 
temperature for the given F (from Fig. 19 or 20 or from in- ' 
terpolation of the models) and has corrected for the 27-day 
variation (formula 12)), one can correct for the geomagnetic 
activity effect by means of formulae (14) or (15) in the 
following way. The correction to the . exospheric temperature is 


AT = aT(A ) [from equ. ;(i4) or (15)]-AT (A [from Table 3]) 
If the temperature increase due to the geomagnetic 

' ■ ■ ■ - - t . 

• ^ K ■ _ . , . V ' • ■ . 

activity effect is to .be known with a higher time resolu- 
tion than one day, one has to use the corresponding three- 

hourly index a instead of A in formula (16). In this 
P P 

case, however, one has' to allow for the reaction time of 

\ 

the upper atmosphere to changes in a . From various storms 

P 


(16 )i 


\ 
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it is found that this reaction tirae i.s approximately 6 

hours. Thus one ought to use the a value for 6 hours 

P 

earlier than the time for which the exospheric temperature 
is to be derived. I'Jhether the delay time depends also on 
local time and on latitude is not yet known. The possible 
dependence on latitude, however, is believed not to be 
important here, since the models apply to the equatorial 
and temperate zones only. 


c ) The semi-annual varia c i o n 


The semi-annual variation in atmospheric densities 


is now observed ovei 


:i;r.e-per iod of 6 years. The 


height interval of the observations ranges from 200 to 
1600 km. The variation seems to occur v/icn considerably 
more regularity than the similar variation in the geo- 
magnetic indices. The main properties of the semi-annual 
variation in the densities and henceforth in tine r.erapera- 
tures are a pronounced minimum in June-vJuly and a somewhat 
less conspicuous minimum in January. The miaximum in 
September -October usually appears to be somev/hat higher 
than the maximum in March-April. The different depths 
of the two maxima led Paet told ( 1 963 ) to separate this 
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effect into a semi-annual and an annual compousu t . There 
is, however, no general agreement as to whether this sapara 
tion also reflects _tv/o different physical processes, it 
is furthermore still controversial v/hether the sarai-annuai 
variation xs related to the solar v/ind as an energy source 
or to a systematic global convection pattern- 

Despite these open 

problems it is, howex'er, fait that the statistical evidence 
is already sufficient to provide a rough correction fonnu la 
again for the exospher j c temperature. It will be vary 
important to prove whether the essential shape of the 
semi-annual variation remains the same for the forthcovaing 
increasing pliase of solar activity in 1965-1960. 

The models presented here have bean adjusted to tha 
semi-annual varie.tion in a v;ay that they represent the 
atmospheric parameter at the maximum of the effect in' 
March-April and an average of tha parameters over tha 
usually higher fall maximum, averaged for the months 
September, October, November. For these reasons t'ne 
correction formula for the semi-annual variation to be 
applied here will differ by an additional quantity from 
those given by Paetnold ( 1964 ) or Jacchia Ov?ftr) • ■'he 
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certain 


difference betv/een the exo.-pheric temperature for a 

date and the temperature calculated frora F, F and A is 

P 


-/^0.39 + 0.15 sin (2n ^--=-^-2, 

355 V 


sin(4TT - 0.30]f (17) 


Here d is the date in days counted from January 1 of each 


year. 


In Fig. 13 


s reviev? (196-) 

the semi-annual variation -.'o- 


^ henavior of 


■ yoci unti.' 


It can b 


seen that the oscillc- 


-inoo IS plotted.' 

r.o;:i exists in each of 
these years in general Iv the se"^ 


seen that 


*-ne time of the maxi'.r.a or ri"- 


ir IS also 


va.rv ov 


^ i"'-' Oc; a) c O': ' ~ d: '* h _ •) 

►-•—-*^ 1 ., 1 iOW 


well over d laonth. Thi 

well this effect can be taken ^nto ^ 

-- ^ .. i ^ ct .. 1 o V / n ci. u 'c a t 

accuracy of the reduction procedure is. 



In order to find the atmospheric model parameters 
(temperature, density, pressure, scale heights, number 
densities or mean molecular weights) for a certain date 
and certain hour of local time for a location in the 
equatorial or temperate zone of the earth at a height be- 
tween 120 km and 800 km, one has to follow the following 
steps: Calculate the average of the solar 10.7 cm flux F 

over a period of 5 solar rotations. (5— monthly averages 

for 1953 - 1963 are given in Table 2. They might con- 
veniently be used.) Then take the solar 10.7 cm daily 
flux value P for the day preceding your date. (The 10.7 
cm flux values are measured by the National Research 
Council, Ottawa, Canady. 'They are published in the monthly 
issues on “Solar Geophysical Data" by the U.S. National 
Bureau of Standards, Central Radio Propagation Laboratory, 

Boulder, Colorado.) Next obtain the geomagnetic index a 

• ,, , . p 

for the time 6 hodrs before your time. These can be de- 
rived from the indices . published monthly in the Journal 
of Geophysical Research by using the following conversion 
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o 

o 

w 

Cr 

1- 

lo 


2- 

2o 

2 + 

3- 

3o 

3 -r 

4- 

4o 4 + 

a =0 
P 

2 

3 

4 

5 

6 

7 

9 

12 

15 

18 

22 

27 32 

K = 5- 
P 

5o 

5 r 

6- 

6o 

6T 

7- 

7o 

7 + 

8- 

8o 

8+ 

9- 9o 

a =39 
P 

48 

56 

57 

80 

94 

111 

132 

154 

17 9 

207 

236 

300 400 

Then proceed 

clS 

fol 

iov/3 ; 










1) Fir 

^ 4— 

S u 

find 

the 

nigh 

. 'w L. 

ime exGsp' 

heri 

c ten 

;per 

m(l) 

acure T. . 

04 

(at 04 hours 

local 

time) 

co:c 

x e 

spondi 

,ng 

to F 

from Figs. 

19 or 


20 or by interpolating the models by raeoMS of Table 3. Reduce to 

I 

A =0 by subcractincr r.he aooropriaite value of AT as given in Table 

p - . . 

2) Correct for the 27-dav vc.i'iaticn. The corrected 


nighttime temperature T^”^' is 


/o' 

s a- J 


/ ^ \ 
\ - y 


'0^ 




L '• 


(IS) 


— p 9 2 

10 “ Vv/m c/s.) 


(The solar 10.7 era flux is ex.pressed in units o 
The response tiraa of the ;-,u..'.osphere is £;pproxir.ms.tely accounted 
for by using the T-valua for one cay prior to the date. 

3) Correct for the gec.r-agnexic - -t- 


corrected nighttime cemparature T„ , 


, \ 
\ - j 

.1 


,(^)_ rr,(2) 

■04 ‘04 


+ 1.0 a + 125'"^ 

p 


exp (“o. 03 a )J [ Kj (r9) 
P 


Account approximately tor the reaction time by tah_rig the 


! 


a value for a time 6 hours prior to your time. 

P 

4) Correct for the semi-annual variation. The cor- 

(4) . 

rected nighttime temperature is 

’’of ” ’'of ^et)" 


Here d is the day within the year (Jan. 1 corresponding to^ 


1 ). 


(4) 

5) Interpolate the models corresponding to in 

order to find the corresponding values of the temperature, 
density, pressure, etc. for the local time for which you 
want to know the atmospheric parameters. 

This procedure can, of course, only be approximately 
correct, since it is essentially based on the presumption 
that the height dependences, of the different heating 
mechanisms which cause the three aforementioned effects 
are nearly the same. .Nevertheless, th® P^®“ 

cedure for finding the atmospheric parameters seems to 
provide reasonably useful results. 

From the discussion in the last two chapters, in 
particular on the accuracies of the model data, it is 
evident that the meaning of our set of reference models 
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cannot ba to provide exact values of the atmospheric 
parameters for each moment and location but is intended 
rather to provide a tool for the intercomparison of forth- 
coming measurements v;hich in general will have been taken 
at different times and locations. These coraparisons in 
turn mignc men lead to an improvement of our knov;ledge 
of the behavior of the upper atmosphere and of the im- 
portant physical processes which determine the atmospheric 
structure and its dependence on time with short-period 
and long-terra variations. Wa are hopeful than the models 
presentea here prove to be helpful in this resoect. 
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FIGURE CAPTIONS 


Pig. 

Fig. I 

Fig. 3 
Fig. A 

Fig. 5 
Fig. 6 
Fig. 7 

Fig. 8; 
Fig. 9: 


1: Temperature of the upper atmosphere as a function 

of height in the range from 120 to 800 km for 4 
selected hours of local time (04, 10, 14 and 20 
^t)urs) during a low level of solar ectivity 
(F = 75). (From Model 2) 

2: Temperature versus height (same as in Pig. 1) for 

mean solar activity^as represented by an average 
solar 10.7 cm flux F = 150 (in units of 10"22 
W/m-^ c/s) . (From Model 5) 

I; Temperature versus height (same as in Fig. 1) for 
very high solar activity (F = 250) , (From Model 9) 

: Logarithms (base 10) of the density in the height 

range from 300 to 800 km for 4 selected hours of 
local time (04, ]^0, 14 and 20 hours) during low 
solar activity (F = 75). (From Model 2) 

: Logarithms of density versus height (same as in 

Fig. 4) for mean solar activity (F = 150). 

(From Model 5) 

: Logarithms of density versus height (same as in 

Fig. 4) for very high solar activity (F = 250). 
(From Model 9) 

: Temperature versus local time for a height of 300 

and of 600 km above the Earth's ^surface for three 
levels__of solar activity: low (P = 75) - Model 2, 

mean (F = 150) - Model 5, and very high (F = 250) 

~ Model 9. 

Log density versus local time for every 100 km 
from 300 to 800 km for low solar activity (P = 75) 
(From Model 2) , 


Log density versus local ^ime (same as in Pig. 8) 
for mean solar activity (F = 150) . (From Modal 5) 



Fig. 10 
Fig. 11: 

Fig. 12; 

Fig. 13: 

Fig. 14: 
Fig. 15: 

Fig. 16; 
Fig. 17: 

Fig. 18: 
Fig. 19: 


Log density versus local time (same as in Fig. s) 

moLI'T 

Mea.-. x-nolecular weight versus height in the 
from 120 to 800 Icm for 4 selected hoL"s oj [oca! 

liveis”"*- ‘>'® three 

^eL ?°! [sor-'Hilt - Model 2 , 

- Model 9. ■ = 250) 

10, 14 and 20 ,ours-r!o1 

2, mean 

i‘iOu.ex 9. 


Logarithms (h 


^2' ^2' 


■? r' \ 

-O; O] 


rnd H veious 


number densities of 

from 120 to SOO fr*- height in the range 

-e .Vi.. 1 . 01 . C-i nours local time and 

low solar accrvity (}? c= 75 ^ fT?-,-,—-. -■> 


Sc^iT.B as Fag. 13 
(Fro.m Modal 2) 


aou.rs local time. 


Same as Fia. 1 3 b-v - 

— ^ y 1- iT is. ' a r* c* /-^ i U5 -T - ^ • 1. 

^oO) . .ina local time is 04 hours fF-or-i 
Model 5) noui.s. (Fx.om 


Same as rig. 15 , 
(From Model 5) 


ror 14 hours local time. 


S^ie as Fig. 13, hub for verv hie 

(r = 2d 0) . The local time is 0-- 
Model 9) 


a a o 1 . ar ac t i vi ty 
p.ours. (From 


Sair.e as Fig. 17 , but for 14 ho)-s ^ 

(Fro.m Model 9) - -o._s .oc. .^me. 


The solar-cycle variation of the exos->be-'b- 
temperature plotted as a function of the le'val 
or soiar activity which is represented bv 
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Fig. 19: 5-monthly averages F of the solar 10.7 cm flux, 

(cont.) Here the values are given for every two hours in 
loc^l time from 04 to 14 hours* 

Fig. 20: Same as Fig. 19. Here the exospheric tempera- 

tures are given for every two hours in local time 
from 16 to 02 hours. 

Fig. 21: Comparison of the models 7 (F = 200) and 2 (F = 75) 

with a few observational data intended to' demon- 
strate an example of the agreement. The__observa- 
tional data are the Bonn model 1961 for F = 200 
(Martin et al. 1962) and the data derived by 
Jac^hia and Slowey for March 1963 corresponding 
to F = 78 and 04 hours local time from the 5 
satellites Injun III, Explorer I, Explorer VIII, 
Vanguard II and Explorer IX. The maximum and 
minimum exospheric temperatures are stated on 
the right. 

Pig. 22: Monthly averages of the solar 10.7 cm flux in 

units of 10“22 from the National Research 

Council, Ottawa, Canada for the years from 1947 
through 1963. The open circles represent the 
yearly averages. The bars indicate the scatter 
of the daily values, represented by the maximum 
and minimum value during the months June and July 
of each year. 
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Model 1 


Properties nf the Upper At.rriOS'Dhere as Functions 
of HeiaVit and Local Time for a Level of 
Extreme Low Solar- ActivltY 

The properties of the upper atmosphere are tabulated 

as a function of local time (for every two hours of the 

day) and as a function of height in the range from 120 km 

to 800 km above the Earth's surface. The quantities listed 

are: temperature [ K], density (g cm ), pressure (dynes cm ), 

scale height (km) , mean molecular weight and the number 
—3 

densities [cm ] of N^/ 0 and He. The last column 

contains the number densities of argon (A) in the range 
from 120, kra to 400 km and the number densities of ''.ydrogen 
(H) in the range from 500 km to 800 km. 

On top of each page the following quantities are 
listed: the local time in hours and the diurnal maximum 

and minimum value of the exospheric temperature for this 
particular model. 

The model is valid for the earth's equatorial and 
temperate zones when the level of solar activity is re- 
presented by an average flux F = 65 x 10 W/m c/s of the 
solar radiation at a wavelength of 10.7 cm. 
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The propex*ties of the upper atmosphere are tabulated 
as a function of local time (for every two hours of the 
day), and as a function of height in the range from 120 Tcm 
to 800 km above the Earth's surface. Tne quantities listed 
are; temperature C°K], density (g cm ), pressure (dynes cm 
scale height (km) , mean molecular weight and the nun^ber 
densities [cm~^3 of N , 0 , 0 and He. The last column 
contains the number densities of argon (A) in the range 
from 120 fan to 400 km and the numZoer densities of hydrogen ' 
(H) in the range from 500 kra to 800 Ion. 

On top of each page the follov/ing quantities are 
listed: the local time in hours and the diurnal maximum 

and minimum value of the exospheric temperature for this 
particular model. 

The model is valid for the earth's equatorial and 

temperate zones when the level of solar activity is re~ 

— ”22 2 

presented by an average flux F = 75 x 10 W/m c/s of the 
solar radiation at a v;avelength of 10.7 cm. 
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Model 3 



Properties of __the .U poar Atmosphere as _.FunctlQiis 
of Height and Local Time for a 
L ov;_ So lar Act.l^far 

The properties of the upper atmosphere are tabulated 
as a function of local time (for every two hours of the 
day) and as a function of height in the range from 120 km 

to 800 km above the Earth's surface. The quantities listed 

■**3 * 

are: temperature [°K], density (g cm ) , pressure (dynes cm 

scale height (km) , mean molecular v;eight and the number 

densities [cm”^] of O 2 / 0 and He. ihe last column .. 

contains the number densities of argon (A) in the range 

from 120 Tcm to 400 km and the number densities of hydrogen ' 

(H) in the range from 500 km to 800 km. 

On top of each page the following quantities are 

listed: the local time in hours and the diurnal maximum 

and minimum value of the exospheric temperature for this 

particular model. 

The model is valid for the earth ’ s equatorial and 

temperate zones when the level of solar activity is re~ 

— -22 2 

presented by an average flux F = 100 x 10 W/m c/s of the 
solar radiation at a wavelength of 10.7 cm. 
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Mode l 4 


■E.^.QP.^JTt’-'ies Qf the Upper Atmosphere as Funct.i.ons 
^_ Haj.ght and Local Time for a Level nf 
Splax Activi.tv Intennediate Between 
Low and Mean Activity 

The properties of the upper atmosphere are tabulated 
as a function of local time (for every two hours of the 
day) and as a function of height in the range from 120 hm 
to 800 km above the Earth's surface. The quantities listed 
are: temperature [°K], density (g cm ^) , pressure (dynes cm 

scale height (km) , mean molecular weight and the numbei' 

-3 

densities [cm ] of N^/ 0^, O and Ke. The last column 
contains the number densities of argon (A) in the range 
from 120 kra to 400 km and the number densities of hydrogen 
(H) in the range from 500 km to 800 km. 

On top of each page the following quantities are 
listed: the local time in hours and the diurnal maxiraum 

and minimum value of the exospheric temperature for this 
particular model. 

The model is valid for the earth’s equatorial and 
temperate zones v/hen the level of solar activity is re~ 
p^esenced by an average flux P = 125 x 10 W/m^ c/s of the 
solar radiation at a wavelength of 10.7 cm. 
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OT-r.r>prtie3 of the Upper Atmosphere •; ..... 

nf Height and Local Time for a Level or 
Maan Solar Act.iv5.tv 

The properties of the upper atmosphere are tabulci.. 

as a function 6f local time (for every two hours of the 

day) and as a function of height in the range from 120 k 

to 800 hm above the Earth's surface. The quantities lis 

are: temperature [°K], density (g crn , pressure (dyne 

scale height (km) , mean molecular weight and the number 
-3 ' ' 

densities [cm ] of K, , 0., 0 and He. The last column 
contains the number densities of argon (A) in the range 
from 120 km to 4-00 km and the number densities of hydrog 
(H) in the range from 500 kra to 800 km. 

On top of each page the following quantities are 
listed: the local tirae in hours and the diurnal maximum 

and minimum value of the exospheric temperature for this 
particular model. 

The model is valid for the earth ' s equatorial and 
temperate zones when the level of solar activity is re- 
presented by an average flux F = 150 x 10 W/m c/s of 
solar radiation at a wavelength of 10.7 cm. 
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Model 6 


Properties of the Upper AtmosT:?here as Functions 
Qf- Height and' Local Time for a Level of 
Solar Activity Intermediate Between 
Mean and Hich Acti^rjtv 

The properties of the upper atmosphere are tabulated 
as a function of local time (for every tv;o hours of the , 
day) and as a function of height in the range from 120 km 

to 800 km above the Earth’s surface. The quantities listed 

”3 — 

are: temperature [ k], density (g cm ), pressure (dynes cm 

scale height (km), mean molecular weight and the number 
-3 

densities [cm ] of N^, 0^. 0 and He. The last column 
contains the nvimber densities of argon (A) in the range 
from 120 km to 400 km and the number densities of hydrogen 
(H) in the range from 500 km to 800 3cm. 

On top of each page the following quantities are 
listed: the local time in hours and the diurnal maximum 

and minimxim value of the exospheric temperature for this 
particular model. 

The model is valid for the earth's eq-uatorial and 

temperate zones when the level of solar activity is re- 

-22 2 

presented by an average flux F = 175 x 10 W/m c/s of the 
solar radiation at a wavelength of 10.7 cm. 
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HXcrh Solar Activity 


The properties of the upper etraosphere ere tabulate 
as a function of local tiraa (for every two hours of the 
day) and as a function of "neigl’t in the range frora 120 k 
to 800 kui above the Earth's surface. The quantities lis 
are: teraperaturs [ K], density (g c.v: , pressure (dyne 

scaie height (krr.) , mean rolecular v/eight and the nur.iber 
densities [cva "^] of N^, O and He. The last column 


comcarns me nuno^ 


in itias of argon (A) in the rang 


from 120 kai to -1-00 km and the nuaber densities of hvdrcc,- 


(K) 


On top of each page the roliowing quantities are 
listed: the local tine in hours c.nd the diurnal maximum 


ana manimum varue or one excsu'.ierac 




ror tnis 


particular m: 


The model is v^alic for the earuh's ecuatoriai anc 


temperate cones whan the raver or sorar actrvrty xs re- 

-- -"2 2 

presented by an average flux i = 250 x 10 Vm c/s or 
solar radiation at a wa.valength of 10.7 cr:.. 
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